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Abstract 

Consider N points randomly distributed along a line segment of unitary length. A walker 
explores this disordered medium moving according to a partially self-avoiding deterministic walk. 
The walker, with memory [i, leaves from the leftmost point and moves, at each discrete time step, 
to the nearest point, which has not been visited in the preceding fi steps. We have obtained 
analytically the probability Pn(^) = (1 — 2~ M ) iV ~ M ~ 1 that all N points are visited in this open 
system, with N ^> fi ^> 1. The expression for Pn(p>) evaluated in the mentioned limit is valid even 
for small N and leads to a transition region centered at \x\ = In N/ In 2 and with width e = e/ In 2. 
For fi < hi — e/2, the walker gets trapped in cycles and does not fully explore the system. For 
fi > fii + e/2 the walker explores the whole system. In both cases the walker presents diffusive 
behavior. Nevertheless, in the intermediate regime fi ~ jii ± e/2, the walker presents anoumalous 
diffusion behavior. Since the intermediate region increases as lniV and its width is constant, a 
sharp transition is obtained for one-dimensional large systems. The walker does not need to have 
full memory of its trajectory to explore the whole system, it suffices to have memory of order jjL\. 
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While random walks in regular or disordered media have been thouroughly explored 
deterministic walks in regular ^ and disordered media 0, 0, 0, (| have been much less 
studied. Here we are concerned with the properties of deterministic walks in random media. 

Given N points (cities) distributed in a (^-dimensional space, a possible question to ask 
is how efficiently these cities can be visited by a walker who follows a simple movimentation 
rule. In the travelling salesman problem, one searches for the shortest closed path, which 
passes once in each city. This problem has been extensively studied. In particular, if the 
coordinates of the cities are distributed following a uniform deviate, results concerning the 
statistics of the shortest paths have been obtained analytically a y, a. To tackle this 
problem, one imperatively needs to know the coordinates of all the cities in advance. Global 
system information must be at the walker's disposal. 

Nevertheless, other situations may be envisaged. For instance, suppose that only local 
information about the neighborhood ranking of the current city is at the walker's disposal. 
In this case, one can think of several deterministic and stochastic strategies to maximize the 
number of visited cities, while trying to minimize the travelled distance. 

Our aim is to study the way a walker explores the medium following the deterministic 
rule of going to the nearest point, which has not been visited in the previous fi discrete time 
steps. We call this partially self-avoiding walk of the deterministic tourist walk [lol . 11 . 12]. 



Each trajectory, produced by this deterministic rule, has an initial transient of length t 
and ends in a cycle of period p. Both transient time and cycle period can be combined in 
the joint distribution S^J(t,p). The fi = deterministic tourist walk is trivial, the walker 
does not move at each time step. The transient and period joint distribution is simply 
S^fditiP) = $t,o8p,i, where Sij is the Kronecker delta. With memory \i = 1, the walker must 
leave the current city at each time step and the transient time and period joint distribution 
has been obtained for N 3> 1 [13]. The \i = 1 walks do not lead to the random medium 
exploration. This is due to very short transient times and the tourist gets trapped in pairs 
of cities, which are mutually nearest neighbors. Interesting phenomena occur when fi > 1 
is considered. In this case, the cycle distribution is no longer peaked at p m in = H + 1, but 
presents a whole spectrum of cycles with period p G + 1, N], with possible power-law 
decay [lO, 11, 12]. These cycles have been used as a clusterization method 14], texture 



analysis in images HQ , primate foraging [4j, |6| and thesaurus analysis 



121. 



It is interesting to point out that, for ID systems, determinism imposes serious restric- 
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tions. For any \i value, cycles with period p G [2/i + 1,2// + 3] are forbidden. Additionally, 
for = 2 all odd periods, but p m in — 3, are forbidden. Also, the heavy tail of the pe- 
riod marginal distribuition sffl (p) = J2t Sjfi (t, p) may lead to often- visited-large-period 
cycles 10[. This allows system exploration even for small memory values (/i <C N). 

In this letter, we consider the deterministic tourist walk along ID random systems. 
Through a simple (though not trivial) derivation, we show the existence of a transition 
in the walker's exploratory behavior at a critical memory fix = In N/ In 2 in a narrow mem- 
ory range of width e = e/ln2. This transition splits the walker's behavior in essentially 
three regimes. Also, we show numerically that the final step distribution of these walks 

nn 

is the estimator of the fractionary and non-linear diffusion equation solution |17l . 1181 ] . For 
(/, < fix — £/2, the walker gets trapped in cycles and for (a > (i\ + e/2, the walker visits 
all the cities. In both cases, the walker presents normal diffusion behavior. Nevertheless, 
for (i ~ (ix ± e/2 the walker superdiffuses indicating that (i\ is an optimum memory for 
maximum exploration and minimum movimentation. 

A random static semi-infinite medium is constructed by uncountable points that are 
randomly and uniformly distributed along a line segment with a mean point density r. The 
distances x between consecutive points follow an exponential probability density function 
(pdf): f(x) = re~ rx , for x > and f(x) = 0, otherwise. 

Consider now the tourist dynamics with a walker who leaves from the city Sx, placed at 
the origin of the line segment. The probability Sj^j(n) for the walker to explore n distinct 
cities can be derived as follows. The first (i+l cities are indeed explored, because the memory 
(A prohibits the walker to turn back. Thus, the distances xx, x 2 , • ■ ■ , may assume any 
value in the interval [0, 00). The following steps are uncertain. The walker may move either 
forward to a new city or backward to an already visited city outside the memory window. In 
analogy to the geometric distribution, it is useful to define the exploration probability q~j as 
the probability for the walker to explore a new city at the j-th uncertain step. Thus, q~x can 
be obtained imposing that the distance x^+x must be less than the sum yx = J2k=i x k- Since 
the variables Xxi %2j ■ ■ ■ 1 -Kfi £ue independent and identically distributed with exponencial 
pdf, yx has a gamma pdf. Hence qx = / °° <ij/ 1 r M j/f _1 e~ rj/1 /r(A i ) Jo^ 1 dx fJi+ ire~ rXf ' +1 = 1 — 2 _At . 

The exploration probability q~i for the second uncertain step is not exactly equal to q\. 
Once the distance x^+x must vary in the interval [0,yi], the variables x 2 , x 3 , . . . , are 
not all independent, and consequently y 2 = J2k=2 x k has not exactly a gamma pdf. However, 
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for \i 1, rarely exceeds yx [this probability is just P(x fl+ x > y x ) = 1 — q\ = 2~ M , 
meaning that a weak correlation is present]. Therefore, one can make an approximation 
assuming that 2/2 follows a gamma pdf and considering q~2 ~ q%. The same argumentation 
can be used for the succeeding steps. 

When the city s n is reached, the walker must turn back, finishing the medium exploration. 
Once q~i is taken for all q, the return probability is p = 1 — q = 2 _At and one has: S^j(n) = 
2^(l-2^)"-^ 1 . Notice that r has been eliminated, indicating that the number of explored 
cities does not depend on the medium density. 

We observe that once the memory \x assures the walker to visit at least fi + 1 cities, its 
convenient to define the number of extra explored cities as n e — n — \i — 1. In this variable, 
Sj^l begins at n e = for all fi value and one has E(n e ) = 2 M — 1, which may be interpreted 
as the characteristic reach of the walk, and Var(n e ) = 2 2fJi — 2 M . 

The finite disordered medium is constructed by N points, whose coordinates are randomly 
(uniform deviate) generated in the interval [0, 1]. The exploration and return probabilities 
obtained for the semi-infinite medium may also be applied to this finite medium. The 
equivalence between these two media can be shown restricting the semi-infinite medium 
length to the first iV points and normalizing it to fit in the interval [0, 1]. For both media, 
the abscissas of the ranked points follow a beta pdf. 

The probability Pjv(a*) f° r the exploration of the whole iV-point medium can be derived 
noticing that the walker must move forward N — (jtt + 1) uncertain steps and, when the last 
city is reached, there is no need to impose a return step. Therefore 



which is plotted in Fig. [TJa). 

In Fig. [H(b) one sees that the probability of full medium exploration increases rapidly from 
to 1, in a well defined transition region. Considering N 3> Eq. [T]may be approximated 
to P n (/j) = (1 - 2-^. To obtain the critical memory fix, one considers the inflexion point, 
leading to 



which is the number of bits to represent the system size. For all N, the curve slope at /ii is 




(1) 



log 2 iV 



(2) 
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In2/e [Fig. 11(b)], the transition region has a constant width 

£ = -^-^3.92. (3) 
ln2 v ; 

This indicates that as N increases, p\ slowly increases, but its deviation is independent of 
N, so that a sharp transition is found for N ^> 1. The comparison between the analytical 
result and computer simulation is depicted in Fig. QJb). Observe that the approximation 
employed leads to satisfactory results even for small N. 

Up to now, we have studied the deterministic walks in ID disordered media with open 
boundary condition. Now, let us focus on the diffusion process. This tells us about the 
medium exploration bulk characteristics. Using periodic boundary conditions for each map, 
the walkers start to move, following the tourist rule, from the most central point of the 
interval [—1/2,1/2]. After iV steps, the walkers enter the periodic part, which can be 
interpreted as the steady state of the exploratory behavior (no new visited cities). 

For a given reduced memory jx = p.jp.\ and after t = N steps, the final position pdf 
p(x,t) has been estimated by P/x(x) = n(x)/(MAx), where M is the number of maps and 
n(x) is the number of travelers with final position between x and x + Ax. The normalized 
distribution is shown in Fig. [21 where one sees the huge dispersion around p = 1. 

We have shown numerically that Pp,(x) ~ p(x,t) is the solution of the fractionary and 



non- linear difusion equation [17. 18]: 



d tP (x,t) =Dd 2 x ^p 2 -\x,t) . (4) 

Here £ is the fractionality parameter, which is related to the step-length pdf second-moment 
divergence, q is the non-linearity parameter, which expresses the influence of correlation 
between step lengths and D is the diffusion coeficient. 

For — fix\ > e/2, p(x,t) has gaussian (£ = 0, q — 1) shapes indicating normal diffusion. 
Further, as \p — pi\ increases, the variance a 2 tends to zero and p{x,t) becomes a Dirac 
delta function: ;1 (p) = 5 Pi2 or S^'^p) w 5 p>N . 

Anomalous diffusion occurs for — p x \ < e/2. For \p — pi \ smaller, but closer to e/2, the 
walks have finite second moment. The q parameter becomes important and the solutions of 



Eq. HI with £ = 0, are the generalized g-gaussians [181 ]: G q (x) = r ) q e q {—x 2 /a 2 q )/^Jd^K, where 
e q (x) = [1 + (1 — q)x] 1 /( 1 ~ qS) is a generalization for the exponential function, a 2 (t) = [2Dt(2 — 
<?)(3 — g)-y4~ 1 ] 2 /( 3 ~'?) is related to the variance [a 2 = a 2 /(5 — 3q)] of p and normalization 



is possible only for q < 3: 7 9 = yj\q - T\{T[l/(q - 1) - 6]/T{(3 - q)/[2(q - 1)] - S}} 1 - 25 , 
where 5 = 1 for q < 1, which expresses sub diffusion and 5 = for 1 < q < 3 expressing 
super diffusion. If q — > 1 + , then 7 — > 1, and one has normal diffusion. The variance is finite 
only for q < 5/3, otherwise the g-gaussians give place to Levy a-stable pdf (£ 7^ 0), which 
are the solutions of Eq. @] around \li — Lii\ ~ 0. 

To fit P/i{x) to the theoretical model p(x, i), we have computed its variance for each ft and 
used the maximum likelihood estimation method to find q (Inset of Fig. [3]). One sees that a 2 
diverges in the transition region presenting a heavy tail if compared to a normal distribution 
(Fig. [2]). For ft < log 2 and ft > 2.2577, Pp,(x) shows a gaussian behavior (q = 1 and £ = 0). 
This is like a random walk model with independent steps and finite variance. In one hand, 
for ft < I, this is due to the short transient length and cycle period. The walker excursion 
is limited to the region near the origin. On the other hand, for ft > 1, this occurs since all 
cycles have a period of N. Thus, after N steps, the walker returns to the starting point and 
the variance vanishes. This implies that the starting point coordinate follows a gaussian pdf. 
In the g-gaussian model, a 2 diverges at q = 5/3. Despite this result being in accordance 
with the estimated q values, the model does not fit the data for 2 log 2 < ft < 2.1072. In 
the region 2 log 2 < ft < 3 log 2 the distributions are Levy a-stable, because their tails are 
heavier than g-gaussian ones. For 3 log 2 < ft < 2.1072, the distributions are almost uniform, 
with a clear central peak. This behavior is probably due the periodic boundary conditions. 

In short, for < ft < log 2 and ft > 2.2577, the tourist follows a traditional diffusion 
process (gaussian solution). For log 2 < ft < 2 log 2, and 2.1072 < ft < 2.2577, it is a non- 
linear superdifusion process (anomalous diffusion with g-gaussian solution). In this region 
the anomalous behavior is obtained since the steps are not all independent. For 2 log 2 < ft < 
2.1072, the process presents another kind of superdiffusion, where the anomalous behavior 
is strongly due to the fact that the variance is not finite. In the region 2 log 2 < ft < 3 log 2 
the distributions are Levy a-stable. 

Finally, the tourist rule can be relaxed to a stochastic walk. Thus, the walker goes to 
nearer cities with greater probabilities. These probabilities are given by an one-parameter 
(inverse of the temperature) exponential distribution. This situation has been studied for 
fi = [3| and fi = 1 {20I and we have detected the existence of a critical temperature 
separating the localized from the extended regimes. It would be interesting to combine 
both stochastic movimentation (driven by a temperature parameter) and memory (//) in 
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the tourist walks to perform full compromise between medium exploration and distance 
travelled. 

In conclusion, if the walker has critical memory /ii the solutions of Eq. H] have infinite 
variance. This favors the walker to explore the whole system with minimum displacement, 
indicating that fi± is an optimum memory for an exploratory strategy. It is intriguing that 
a simple deterministic system as this one can present a complex behavior given by the 
complete fractionary non-linear difusion equation and that a small memory value allows the 
global medium exploration. 

The authors thank N. A. Alves and F. M. Ramos for fruitful discussions. ASM acknowl- 
edges CNPq (305527/2004-5) and FAPESP (2005/02408-0) for support. 



* Electronic address: cesartercariol@gmail.com 
' Electronic address: caminhos'rsg@yahoo. com.br 



■f Electronic address: | asmartinez@ffclrp.usp.br 

[1 



[2] 
[3] 

[4 

[5 
[6 

[7] 
[8] 

[9 
[10 

[11 
[12 



D. S. Fisher, Phys. Rev. A 30, 960 (1984). 
H. Freund and P. Grassberger, Physica A 190, 218 (1992). 
L. A. Bunimovich, Physica D 187, 20 (2004). 

D. Boyer, O. Miramontes, G. Ramos- Fernandez, J. L. Mateos, and G. Cocho, Physica A 342, 
329 (2004). 

D. Boyer and H. Larralde, Complexity 10, 52 (2005). 

D. Boyer, G. Ramos-Fernandez, O. M. J. L. Mateos, G. Cocho, H. Larralde, H. Ramos, and 

F. Rojas, Proceedings of the Royal Society B-Biological Sciences 273, 1743 (2006). 
A. G. Percus and O. C. Martin, Phys. Rev. Lett. 76, 1188 (1996). 

N. J. Cerf, J. H. B. de Monvel, O. Bohigas, O. C. Martin, and A. G. Percus, J. Phys. I (France) 
7, 117 (1997). 

A. G. Percus and O. C. Martin, J. Stat. Phys. 94, 739 (1999). 

G. F. Lima, A. S. Martinez, and O. Kinouchi, Phys. Rev. Lett. 87, 010603 (2001). 

H. E. Stanley and S. V. Buldyrev, Nature (London) 413, 373 (2001). 

O. Kinouchi, A. S. Martinez, G. F. Lima, G. M. Lourengo, and S. Risau-Gusman, Physica A 
315, 665 (2002). 



8 



[13] C. A. S. Tercariol and A. S. Martinez, Phys. Rev. E 72, 021103 (2005). 

[14] M. G. Campiteli, P. B. Diniz, O. Kinouchi, and A. S. Martinez, Phys. Rev. E 74, 026703 
(2006). 

[15] A. R. Backes, O. M. Bruno, M. G. Campiteli, and A. S. Martinez, Lect. Note Comput. Sci. 
4225, 784 (2006). 

[16] M. G. Campiteli, A. S. Martinez, and O. M. Bruno, Lect. Note Comput. Sci. 4140, 159 (2006). 
[17] R. Metzler and J. Klafter, Phys. Rep. 339, 1 (2000). 
[18] C. Anteneodo, Physica A 358, 289 (2005). 

[19] A. S. Martinez, O. Kinouchi, and S. Risau-Gusman, Phys. Rev. E 69, 017101 (2004). 
[20] S. Risau-Gusman, A. S. Martinez, and O. Kinouchi, Phys. Rev. E 68, 016104 (2003). 



9 



(a) 
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FIG. 1: (a) Probability of full exploration (Eq. [1]) as a function of \x and N for ID deterministic 
tourist walks. One sees the abrupt transiton between the localized (bottom) and extended (upper) 
regimes, (b) Projections of Eq. [T]for some fixed N values. Analytical results are satisfactory, when 
compared to numerical simulation (M = 10 000 maps for each N and fi values), even for small N 
and /i values. Error bars are smaller than symbol size. The transition points \i\ are given by Eq. [21 
which are weakly dependent on iV but all of them have the same constant dispersion e ~ 4 (Eq. [3]). 
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FIG. 2: Distribution of the final position after t = N steps in the ID tourist walk as a function 
of Numerical simulations have been perfomed with Ax = 1/100, M = 10 5 maps with N = 10 4 
cities each using periodic boundary conditions. The memory \i varied from 1 to 50. 
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FIG. 3: On the lefthand side axis one reads the variance a 2 (•) of curves Pn(x) (estimator of p) as 
function of normalized memory jj,. Observe that a 2 diverges near the critical memory. This region 
is equivalent to the smallest amplitude, but with maximum dispersion of P^(x), in Fig. El On the 
righthand side there are the values of q (A) as function of \x. Notice q ~ 5/3 in the region where 
a 2 diverges. Inset: Pu(x) as a function of the final walker positon for some memory values. The 
displayed curves are g-gaussians defined in the text. 
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